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In the skutterudite compounds the anharmonic 'rattling' oscillations of 4f-host ions in the sur- 
rounding Sbi2 cages are found to have significant influence on the low temperature properties. 
Recently specific heat analysis of Pr(Osi_ a: Ru a; )4Sbi2 has shown that the energy of crystalline elec- 
tric field (CEF) singlet-triplet excitations increases strongly with Ru-concentration x and crosses the 
almost constant rattling mode frequency ujo at about x ~ 0.65. Due to magnetoelastic interactions 
this may entail prominent nonadiabatic effects in inelastic neutron scattering (INS) intensity and 
quadrupolar susceptibility. Furthermore the Ru- concentration dependence of the superconducting 
T c , notably the minimum at intermediate x is explained as a crossover effect from pairforming 
aspherical Coulomb scattering to pairbreaking exchange scattering. 
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I. INTRODUCTION 



The tetrahedral rare earth skutterudite compounds RM4X12 (R= rare earth, M = Ru, Fe, Os and X = P, As , Sb) 
are found to exhibit a great variety of electronic ground states driven by 4f electron correlations jl| . Metallic heavy 
electron and mixed valent as well as Kondo semiconducting behaviour may be found. The 4f ground state may show 
quadrupolar or more exotic multipolar order 043- ^ n PrOs4Sbi2 and substituted compounds superconductivity 
appears 0-0 which is possibly of strongly anisotropic nature • The latter has recently also been found in 
PrPt4Gei2 [Hi ■ It has been proposed [TJ [l2j that low-energy crystalline electric field (CEF) excitations play an 
important role in the formation of Cooper pairs and likewise in the quasiparticle mass enhancement (T3 . Ihlj . These 



only slightly dispersive excitations were identified in INS experiments [Iff in PrOs4Sbi2. Furthermore NMR/NQR 
flfjf and ultrasonic [13, [3 measurements give evidence for the importance of anharmonic 'rattling' oscillations of 4f 
host ions in the cages formed by surrounding Sbi2 icosahedrons on the low temperature properties. 

Recent systematic specific heat measurements in the substitution series Pr(Osi_a;Ru x )4Sbi2 have shown [2(| that 
both CEF excitations and local rattling mode give contributions in addition to the usual Debye part. From the 
analysis Miyazaki et al were able to obtain the dependence of the low energy singlet-triplet CEF excitation A(x) and 
the rattling phonon frequency of the Pr host uio(x) on the Ru - concentration x. It was found that A(x) increases with 
Ru content (x) from A(0) ~ 8 K to A(l) ~ 84 K. For x c ~ 0.65 the triplet excitation energy crosses the oscillator 
energy ujq ~ 45 K of the low energy anharmonic rattling phonon associated with the Pr oscillations in the cage. 
The latter is almost independent of the Ru concentration x. Earlier results @ showed that the critical temperature 
T c (x) varies between T c (0) = 1.85 K and T c (l) = 1.20 K and exhibits a minimum in between with T c (x c ) = 0.7 
K around the same Ru content x c . This raises the question whether there is a connection between the crossing ef- 
fect and the observed T c minimum and whether the nonadiabatic effects play a role in the appearance of the minimum. 

In this work we discuss a model for the nonadiabatic effects between rattling phonons and CEF singlet-triplet 
excitations in Pr(Osi_ x Ru x )4Sbi2. These effects should be primarily observable close to x c where an anticrossing and 
formation of 'vibronic' or mixed modes is expected. The latter have partly a phononic and partly a CEF excitation 
nature. Although similar effects are known for other rare earth compounds [19( they have sofar not been identified 
in the skutterudite family. The observed A(x), u)o(x) behaviour in Pr(Osi_ 2; Ru 2 ;)4Sbi2 however suggests their 
presence. The mixed mode formation should lead to a direct clear signature in the dipolar INS cross section, magnetic 
and quadrupolar susceptibility as well as rattling phonon spectral function. Indirectly other physical quantities like 
NMR relaxation, resistivity and T c suppression or enhancement should also be influenced. Some of these effects 
will be discussed in the present work within an exactly solvable bosonic model for the two types of coupled excitations. 

In Sec. |TI] the vibronic model is introduced and solved in Sec. Mil The dynamical susceptibilities and associated 
structure functions, in particular the dipolar one relevant for INS are derived in Sec. IIV1 Furthermore we investigate 
the influence of vibronic excitations on the formation or breaking of Cooper pairs and the resulting T c (x) variation 
with Ru content in Sec.[V] Finally Sec. I VII gives the summary and conclusion. 
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II. MODEL DEFINITION 



First we consider the local 4f electronic part. The effect of a tetrahedral CEF on the Pr 3+ has been studied 



in detail in Ref. l2ll . Its main consequence is a mixing of the cubic (O^) T 4 and T 5 triplets to tetrahedral T 



(1,2) 
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triplets which have both dipolar and quadrupolar transitions from the ground state singlet Ti. In the following we 
restrict to Ti and lowest T\ = T t triplet states as shown by Shiina [22j. The singlet-triplet CEF Hamiltonian in 
Pr(Osi_ K Ru x )4Sbi2 may then be written in bosonic form as 

H st = A(x)^2(al l a n + -) (n = x,y,z) (1) 

n 

where A(a;) is the singlet-triplet CEF splitting of P^Osi-^Ru-E^Sb;^ that varies from A(0) = 8 K for the Os 
compound to A(l) = 84 K for the Ru system. The triplet states are linear combinations of cubic Ts and states 
[2lj described by 

|i» = (i-d 2 )2|r 5 n) + d|r 4 n) (2) 

where d characterizes the strength of the tetrahedral CEF part (Appendix [A"f . When the latter is small (d <C 1) then 
r t is close to the nonmagnetic cubic T$ triplet. When the tetrahedral CEF part dominates we haved 2 -» \ (Eq.JMJ) 
and r t is an equal-amplitude mixture of nonmagnetic and magnetic T4 cubic triplets. The |r t n) states are created 
by the bosonic operators a„ (n=l-3) according to |r f n) = ajjr s ) from the singlet ground state [22, [23|. They are 
related to a n (n = x,y,z) through a\ = — [1/ \/2)(a x — ia y ); a-i = a z ; 03 = (l/y2)(a x + ia y ). The singlet-triplet 
system has quadrupolar and , because of the tetrahedral mixing amplitude d, also dipolar matrix elements for inelastic 
transitions. In bosonic representation the dipole operators are given by 

J = b D (a x + al,a y + a y ,a z + a\) (3) 

where br> = 2^/|c? is the dipolar matrix element ~ d. The Ts-type quadrupolar operators O n (n — yz, zx,xy) are 

generally given in terms of the total angular momentum components J n (n = x, y, z). In bosonic representation one 
has 

O yz = JyJ z + JzJy = ibQ(a x - a£) 

O xz = J X J Z + J Z J X = ib Q (a y - a\) (4) 
O xy = J x J y + JyJ x ^ ibq{a z - a\) 



where 6q = -^ 2 =-\/35(l — d 2 ) is the quadrupolar singlet-triplet matrix element which is maximal for d=0, contrary to 
bu- These are the order parameters in the field-induced antiferroquadrupolar phase of PrOs4Sbi2 [HI, [24l - [27j and 
their dynamics corresponds to the excitation of quadrupolar excitons (l5l l23l | . 

Now we introduce the rattling phonon part which may be seen as a low frequency optical phonon corresponding 
to the anharmonic movement of the heavy Pr ion in the wide cage formed by the Sbi2 icosahedron. Such almost 
dispersionless rare earth host modes lying within the acoustic phonon bands are reported in Ref. [28| for the Ce 
skutterudite. They belong to Ti representation of Th and therefore are triply degenerate, corresponding to the three 
Cartesian directions of rattling motion in the cage. Due to the anharmonic potential of the cage the effective rattling 
frequency w e (T) may be temperature dependent, similar as in the ft- pyrochlore superconductor KOS2O6 [29l.l30j. On 
the other hand the low temperature effective rattling frequency wo = uj e (T = 0) is almost independent of Ru content 
x with ujq(x) — 45 K. Therefore, as observed in Ref. [2CJ the singlet- triplet energy A(x) crosses the rattling frequency 
around x c ~ 0.65 (dashed lines in Fig. [1} which is, incidentally, close to the Ru concentration where T c (x) shows its 
minimum. In the quasiharmonic approximation [30( the rattling phonon part at low temperature is given by 

ffr = Wo5>K + ^) (5) 

n 

Where n=l-3 denotes one of the triply degenerate host modes which are created by the phonon operators b^ n . 

The coupling of rattling modes and local CEF excitations (all dispersive effects in phonons and CEF excitations 
are neglected) may be written in terms of displacements and 4f quadrupoles as (3lj 

Hr-if = ,9oQo5>n(*) + bt(i))O n (i) (6) 
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FIG. 1: Mode frequencies as function of Ru content x. At zero coupling 7 (dashed lines) fiat rattling mode and increasing 
singlet-triplet CEF excitation cross around x c ~ 0.65. A finite coupling 7 leads to anti-crossing vibronic modes fi^, (upper 
mode) and f)^, (lower mode) around x c . 



where Qq — (2MNu>o) 2 (M=mass of Pr and N=number of sites i) and go is the coupling constant. Expressing the 
quadrupole operators with singlet-triplet boson operators the total Hamiltonian for each site is given by 

T Wo(&],&n + b n bl) + -A(a+a„ + a n a\) + igob Q (b n a n - bt a t + fct a?j _ b n a\)] (7) 



H = 



This is the bosonic model Hamiltonian used in the following analysis. The first part describes three degenerate rattling 
phonon modes the second singlet-triplet CEF exitations and the last one their local magnetoelastic interactions. 

It is bilinear in the singlet-triplet CEF (a) and phononic (b) boson operators and may therefore be diagonalised. 
For that purpose we write it in matrix form as 
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Here we defined go = go(2Mujo) 2 and go — 6qgo- 



(8) 



III. VIBRONIC EXCITATIONS 



The eigenstates of the model Hamiltonian in Eq. © are the local vibronic modes of Pr 3+ , i.e., mixed rattling 
phonon and singlet-triplet CEF excitations. The mode mixing becomes strong close to the crossing of A(a;) with 
wo at x c - The formation of vibronic modes therefore should influence physical properties, in particular close x c - To 
calculate the vibronic modes we express the model Hamiltonian in 2 x 2 block form: 



H = 



n—x.y.z 



(at a„ ) 



Din D2n 
Dzn Din 



(9) 



where a\ L = (a n b n ) and D 3n = D\ n . The D-matrices are the 2x2 blocks in Eq 



(|8]h This quadratic form can be 

diagonalized by a generalised Bogoliubov or paraunitary transformation 9 [32| with the property 916^ = I where 
/ = diag(ll — 1 — 1). In the following the mode degeneracy index (n = x,y,z) will be suppressed. Applying 
the paraunitary transformation we get D9~ 1 = I9^ 1 L where the column vectors w p (p = 1 — 4) of 9~ x are the 
eigenvectors of the equation (ID — X p l)w p = corresponding to the eigenvalue X p determined by the secular equation 
det(D — X p I) — 0. Furthermore L — diag(X+, A_, — A + , — A_). The vibronic eigenvalues can be obtained as 



A± = 


[e ±-B]' 
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FIG. 2: Contour plots of dipolar Sd(u) (left) and rattling phonon S r (w) (right) spectral functions. For better visibility of the 
vibronic effect we choose the coupling constant 7 = 0.3 and a line width T s /loo = 0.15. 



The (transposed) eigenvector column w\ of 9 1 corresponding to eigenvalue A (A p , p = 1 — 4) is given by 
with a normalisation constant 



(# + A)(^f-A)' # + A' (f + A)0f + A) 



(# + A) 2 (^ 



A) 2 (^f 



A) 2 



4#|A|[(f + A) 2 (^-A) 2 + 4g 2 ff] 



(11) 



(12) 



The paraunitary transformation defines the vibronic normal mode coordinates 7 via T] = 8 1 -f where r\ — (aJ,a) = 
(aJb^ab) are the original boson coordinates. Like the latter, the vibronic normal modes fulfill the bosonic commutation 



relations 
leading to 



= 1 



pp' 



Explicitly we write 7^ = (<fv,i/j', <fi, ijj). They diagonalise the Hamiltonian in Eqs. (|8I9[) finally 



H 



(13) 



2A + and 



Where the triply degenerate (n=x,y,z) normal mode frequencies f2" = f2 s (s = 0, tp) are given by fi^ 
Q^p = 2A_. The relation between the non-interacting ^-bosons and the new normal mode coordinates 7 is explicitly 
given by 



/ a n 

bn 



(■ 



J 




(14) 



where the matrix elements 



\w Xi w X2 w X3 WXi 
are obtained from Eqs. (Ill 1 12[) using the eigenvalues in Eq. (fpQ|) . 



IV. DYNAMICAL SUSCEPTIBILITIES AND STRUCTURE FUNCTIONS 



With the above closed solution all interesting susceptibilities and dynamical structure functions of the model may 
be calculated analytically. The dipolar structure function is proportional to the INS cross section and is therefore the 
most direct means to observe the vibronic modes. Furthermore they may influence the NMR rate which is obtained 
from the dynamical dipolar susceptibility. The dynamics of quadrupole moments and the phonon spectral functions 
may be obtained in a similar way. 
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A. Dipolar susceptibility and INS spectral function 

Using the representation in Eq. ([3]) for the dipole operator the dipolar susceptibility may be expressed as 

/>oo 

XD(u) = ib 2 D dte wt {[(a + <J) u {a + a)U)6 H {t) (15) 



where Oh (t) is the Heaviside function. Applying the paraunitary transformation we obtain 

a + a f = u D ^((/) + t ) + Uof(ip + V' t ) (16) 

with uds (s — <f>,tp) given by 

where we defined 7 = 4 go = 4bQgo(2MuJo)~^ . This leads to the dynamical dipolar susceptibility 

X dH - -blY,^D S Ds^) (18) 

S 

where D s (uj) is the retarded Greens function of normal mode bosons according to 

/- on 
dte^{U + </>% (0 + tf) ])6 H (t) = ^.J_ o2 (19) 
-00 l w 1 20 J "0 

and a similar equation for D^(lu). The corresponding spectral function is given by 

Using the fluctuation dissipation theorem w) = e~" u S s (uj) ((3 = 1/kT) the total dipolar spectral function 
corresponding to Eq. (TT5)) is then obtained as 

S D {u) = 2nbl ^ J^T^0^^ [(«- + - «-) + + (21) 

Here n s = (e^ 3 — 1) _1 is the Bose distribution function. When we include a constant finite linewidth T s for the <fi, ip 
bosons the delta functions have to be replaced by Lorentzians. Then at zero temperature we obtain 

S»{ U ) = HEi M ^ F ^ Awo ' ( W _n')2 + r 2 ^ 

Finally, using D' s (0) — — (2/r2 s ) we obtain the zero-temperature static dipolar susceptibility as 

XD{UJ -° ] - {-A) 4- K-^)2 + 7 2Ac = "AT (23) 

Therefore the static dipolar susceptibility of the vibronic system for 7 ^ is unchanged from the dipolar van-Vleck 
susceptibility of the uncoupled (7 = 0) singlet-triplet CEF states because the formation of vibronic modes involves 
only the quadrupolar CEF excitations. The result in Eq. (f2"2"j) together with Eq. (fTU|) gives the frequency and a;, 7 
dependence of the dipolar spectral function which may be compared with INS results. This will be further discussed 
in Sec. 1223 



B. Quadrupolar susceptibility and rattling phonon spectral function 

The dynamical quadrupolar susceptibility, the phonon Green's function and their associated spectral function may 
be obtained in a completely analogous way. Using the bosonic representation of quadrupolar operators in Eq. ([5]) we 
have 

/>oo 

XQ{u) = -ib 2 Q dte iut ([(a-a^ t ,(a-a^o])e H (t) (24) 
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FIG. 3: The top panel shows the calculated T=0 dipolar susceptibility (which is independent of vibronic coupling 7) in 
comparison to experimental values from Ref. 6. Bottom panel shows the (7 -dependent) T=0 quadrupolar susceptibility. It is 
proportional to the T=0 elastic constant reduction AC44. 



Again replacing the a- bosons with <f>, ip bosons by the paraunitary transformation we get 

a -a) = - <f>') + uq^(V> - ■0 t ) (25) 

where now we have a slightly different 



A (^-» 2 ) 2 



Similar as before the quadrupolar spectral function is obtained as 

,2 02^12 



S Q (u) = 2nbl £ A _ [ K + 1)Siu n s ) + n,f (a, + fi.)] (27) 

and for finite boson line width and in the limit T=0 we likewise obtain 

Sq( W ) = 2Aj ^ . (w _^ )2+r , (28) 

Furthermore the static zero-temperature quadrupolar susceptibility may be obtained as 

XQ(0) - (^J (wo2 _ fi 2 )2+7 2 Awo (29) 

The prefactor is the quadrupolar van-Vleck susceptibility of the uncoupled singlet-triplet states. The quantity in 
Eq. (I29|) depends on 7 and therefore on the mode splitting, contrary to xd(0). It is in principle accessible in 
ultrasonic experiments where it determines the velocity or elastic constant change for T — > 0. The latter is given by 
Ac44(T — > 0) = —344X0 where (744 is the magnetoelastic coupling constant of the C44 transverse mode propagating 
along (001) direction [3l|. 

Now we consider the retarded propagator of the rattling phonon which is defined by 

/>oo 

D r {w) = -i dte wt {[{b + tf) tl {b + b%])6 H (t) (30) 



Applying the paraunitary transformation we may express 

b + tf = u^-cb^+u^-^) (31) 
where u s = i\u s \ and the modulus is now given by 



o w 7 2 Aw 



H ~ n s (^ -n^ + ^Au, (32) 
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This leads to a spectral function of the rattling phonon propagator 

7 2 Aw 



s r (uj) — 2n 



CJ 



(n s + l)5(u) - fig) + n s 5{ui + fi s ) 



n s ( w 2-r>2)2 + 7 2 Awo 

Including the finite line width T s for the normal modes we obtain the zero temperature limit 

A 7 2 Aw r s /vr 



S r (u) = 2ir 



n s k 2 -o 2 ) 2 + 7 2 a Wo (w 



n. 



r 2 



(33) 



(34) 



The rattling phonon spectral function in Eq. (1341) is complementary to the dipolar and quadrupolar spectral function. 
In our localized model they are momentum independent. However in the INS cross section the latter is multiplied by 
the square of the electronic form factor F(Q) of the 4f shell which decreases with |Q| while the former is multiplied 
by |Q| 2 . Therefore the dipolar excitation may be seen at small and the rattling phonon part at large total momentum 
transfer. 



C. Numerical results for spectral function and static susceptibilities 

The basic feature of the vibronic mode formation is shown in Fig. [T] At the crossing of the bare (7=0, dashed 
lines) rattling mode u>a(x) and Aq(x) a repulsion takes place for finite magnetoelastic coupling 7 and anti-crossing 
mixed modes (full lines) Q$ (upper mode) and (lower mode) are formed. Their splitting increases with coupling 
strength. At the crossing where luq(x c ) — A(x c ) we have 6 — Cl^ — Qw, — 7 as long as 7/wo "C 1. For 7 = 0.25 the 
splitting is still moderate enough to be compatible with the mode energies determined from specific heat analysis 
|20j. The determination of the mode splitting and f2 s (x) (s = 4>,4>) requires the investigation of the dynamical 
magnetic and phononic structure function in INS experiments. The former should be proportional to Sd{^>) and the 
latter to S r (ui). These functions are calculated from Eqs. (|22l34j) . respectively, and are shown in Fig. Away from 
the anti-crossing region jSd(w) has appreciable intensity only around the bare CEF excitation A(x) and S r (u) only 
around the bare rattling phonon frequency ljq. In the anti-crossing region Sd,t(^) have equal intensity at both split 
modes Q s (x). Observation of this feature by future INS experiments would directly confirm the vibronic mixed mode 
formation in Pr(Osi_ 2: Ru a ;)4Sbi2. 

The static dipolar and quadrupolar susceptibilities Xo,r( w ) are a ls° accessible in experiments. However the former 
does not show an effect of the mode coupling 7 but remains the unrenormalised singlet-triplet van Vleck susceptibility 
(Eq. (I23|l ). Therefore no information on the mode splitting can be gained from it. The behaviour of Xd(%) (T — ¥ 
0), using a constant matrix element bp in comparison to experimental data from Ref. @ is shown in the top Fig. [3] 
Up to x ~ 0.4 the behaviour is in agreement however for larger x the experimental values show no further decrease. 
This may be in part due to the increasing relative importance of higher levels when the singlet-triplet A(x) increases 
(Fig.[IJ. Furthermore if the tetrahedral CEF part increases with x the dipolar matrix element bo ~ d will also increase 
leading to a larger \d at higher x. There is evidence from the superconducting pair breaking behaviour discussed in 
the next section that this is indeed the case. The quadrupolar (Fs-type) susceptibility may be obained from the T=0 
suppression of the appropriate (C44) symmetry elastic constant. The suppression gets larger for increasing vibronic 
coupling 7 at small x. Therefore measuring Ac44(T — ¥ 0;x) may be used as an indirect means to determine the 
coupling strength. 



V. SUPERCONDUCTING PAIR FORMING AND PAIR BREAKING BY VIBRONIC EXCITATIONS 

The effective pairing interaction for the formation of Cooper pairs in skutterudites consists of three contributions: 
(i) harmonic phonons, (ii) local CEF excitations, (iii) low energy rattling (anharmonic) phonons. The former two 
have been included in the model for Lai„ a Pr y Os4Sbi2 [HI (y = Pr-concentration) and the latter were proposed 
for the pyrochlore superconductor KOS2O6 in Ref. H^. However the NMR relaxation [lj|, ultrasonic experiments 
(l7| and specific heat measurements [2(| have indicated that rattling phonons are also present in La, Pr-skutterudite 
compounds and therefore may contribute to the effective pairing mechanism. In fact the INS experiments in 
CeRu4Sbi2 have shown Q the existence of a flat optical phonon mode at wo = 6 meV within the acoustic phonon 
band which shows little hybridisation with the latter. The flat optical mode was interpreted as a (Ce-) guest (rattling) 
mode within the rigid cage structure formed by Sbi2 icosahedrons. This mode is anharmonic and the effective 
frequency u> e (T) is temperature dependent. Around T c however this may be neglected since io e (T) has reached its 
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FIG. 4: Crossover from pair forming (T c > T°) at x =0 with pq > po to pair breaking at x=l with pq < pd (T c < T^). At 
x=0 pd/ Pq = 1/7 as in Ref. [ll| which produces the experimental enhancement T c /T ( ?=2.5. At x=l pd/pq = 1-43 is adjusted 



to the experimental depression T c /T,?=0.34. 



low temperature asymptotic value u>o- In the case of Lai_ y Pr a Os4Sbi2 (x=0) which was studied in Ref. |ll| one has 
A <C uq, therefore the interaction of rattling phonon and singlet-triplet excitation may be neglected. This may not 
hold in Pr(Osi_ a: Ru x )4Sbi2 for general x because A(cc) crosses luq at x c ~ 0.65. Therefore the nonadiabatic vibronic 
spectral function should be used in modeling the effective pairing interaction for arbitrary Ru content x. We mention 
that 'nonadiabatic' refers to the localised 4f electron- phonon interaction, not to the conduction electron- phonon term. 



Experimentally it was found early in Ref. that T c (x) has a minimum close to x c ~ 0.65. Therefore the question 
arises whether this is tied to a suggested vibronic mode splitting or to a different origin. The most convenient starting 
point for a theoretical description is the T°(x) background variation in La(Osi_ K Ru 2: )4Sbi2 which is determined 
by the harmonic and rattling phonon mechanism. The microscopic model behind will not be further specified. The 
symmetry of the superconducting order parameter in the skutterudites is presumably of (anisotropic) extended s-wave 
tyP e [S El- F° r the present purpose we ignore the superconducting gap anisotropy. The presence of 4f states in 
Pr(Osi_ x Ru x )4Sbi2 leads to a scattering of conduction electrons from singlet-triplet CEF excitations. This modifies 
the pair amplitude and changes the background T°(x) of the La compound to a renormalised T c (x). This process is 
due to exchange and aspherical Coulomb scattering of conduction electrons from the 4f shell given by 



H sf = -lac Y /«( ( l) q« c L c k+q ( T - (gj 



17 n ^q« c k(T C kH 



(35) 



qk,ner 



qk,ncr 



Here O qn = ^2 O n (i) exp(iqRi) (n=yz,zx,xy) and J qn = J2q Jn(i) exp(iqRi) (n=x,y,z) are quadrupolar and dipolar 
operators, respectively. Furthermore gj — 4/5 is the Lande factor and / nq = q y q z ,<lz<lx,<lx<ly are quadrupolar form 
factors (q = q/|q|). The principal effect of H s f on superconducting properties of 4f compounds with CEF splitting 
has been investigated in Ref. l33l It was found that for singlet superconductors aspherical Coulomb (quadrupolar) 
scattering which supports pair formation and enhances T° because O n (Eq. is even under time reversal. In 
contrast the exchange term leads to pair breaking and reduces the background T° because J n is odd under time 
reversal. For the case of having only a twofold Kramers degenerate ground state level the latter is described by the 
well known Abrikosov-Gorkov 3J] theory. The modified T c (x) of Pr(Osi_ a; Ru :r )4Sbi2 includes both effects because 
the singlet-triplet excitations have dipolar as well as quadrupolar character due to the tetrahedral CEF (Sec. [IT]). I n the 
present case their magnetoelastic interaction with rattling phonons leads to vibronic excitation modes with modified 
dipolar and quadrupolar matrix elements. Generalization of the expressions for pure CEF systems in Refs. {ll|, IHF 
to vibronic excitations leads to an equation for the renormalised T c given by 



8,Tc_ 
ttVTP 



In 



pq 



J2 u 2 Q MF(^r) + PD J2 



2T C 



(36) 



Here u 2 Ds ,Uq s are given by Eqs. (|17I26[) . The dimensionless vibronic pair forming and breaking strengths Pq{x) and 
Pd{x) for the quadupolar and dipolar conduction electron scattering channels due to H s f are given by 



Pq( x ) = VPQi x ) 



27TN F P ac (P) 

T c 0(0) 



36 2 Q (0) 
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FIG. 5: Dependence of critical temperature on x (Ru-content). Top: Dash-dotted line shows interpolated T°(x) curve for 
La(Osi_ a; Ru a ;)4Sbi2. Only T°(0) = 0.74 K and T°(l) = 3.58 K (full circles) are experimentally known values. Full line shows 
the calculated renormalisation factor T c (x)/T°(x) (with 7 = 0). Note that T°(0) < T c (0) = 1.85 K but T°(l) > T c (l) = 1.20 
K. Bottom: T c (a;) for Pr(Osi_a;Rua;)4Sbi2. Full line is obtained from top panel. The experimental T c data (squares) are taken 
from susceptibility results in Ref. [g. The minimum is due to a decreasing T c (x)/Tc(x) and an increasing T°(x). Dashed line 
(7 = 0.15) shows that effect of vibronic level splitting on T c (x) is small. 



Pn(x) = VPnW-^^y Pn(x) = ^ 3^(0) (37) 



where y is the Pr content. In the present case of Pr(Osi_ 2; Ru a ;)4Sbi2 we have y = 1. Furthermore is the conduction 
electron DOS of La(Os 1 _ ;E Ru a ;)4Sb 1 2 at the Fermi level and t a c (x) = T ( ?(x)/T ( ?(0) the normalised background transition 
temperature of La(Osi_^Ru a: )4Sbi2 (y = 0) with T c (0) — 0.74 K. The quadrupolar and dipolar matrix elements 6q, &d 
are given in Appendix 1X1 and (f 2 ) is the Fermi surface averages of the quadrupolar form factors (independent of n). 
The functions F, G in Eq. (|36[) which correspond to pair formation and pair breaking respectively in principle depend 
explicitly on temperature via the thermal occupation of excited vibronic levels [33| . However in the present case we 
have A(x)/2T®(x) ^ 1 and therefore £l s (x)/2T®(x) ^> 1 for all Ru concentrations x. Then we may use the low 
temperature limit for F, G in which case we have, defining x s — f2 s /2T c : 

F(x s ) = -L + Sl (x s ) 

Xs 

G(x s ) = —-S 1 {x s ) + S 2 (x s ) = -F{x s ) + S 2 (x s ) (38) 

Xs 

Here Si^(x s ) are combinations of digamma functions derived in Ref. |33| and given in Appendix [Bl for completeness. 
The dimensionless pq(x), pd(x) parameters in Eqs. (|36l37p depend on the Ru content x via three factors: (i) the 
quadrupolar and dipolar matrix elements bi(x) (i = Q, D) which are determined by the x-dependent CEF parameters 
xcf{%) and ycF(x) (Appendix |A"|) . (ii) the dimensionless interaction constants f>i(x) depend on x through conduction 
electron DOS Nf and possibly also through the interaction strengths l aCl l ex . (iii) the background normalised 
transition temperature t° c (x). 

A model for these Ru concentration dependences is needed as input to calculate the T c (x) for the Pr compound 
from Eq. ((3!)]) or T^(x)/T^(0) normalised to the La(Osi_ x Ru x )4Sbi2 value. The model for bi(x) (i = Q,D) and 
Pi(x) is described in Appendix |A"1 The background Tc( x ) i s determined experimentally only for the stoichiometric 
cases (x — 0, 1). The interpolation in Fig. [5] is used for intermediate concentrations. The resulting x-dependence of 
pair- forming and breaking parameters Pq{x) and Pd{x) obtained from Eq. (|37[) is shown in Fig. U and fit procedure 
and parameters are explained in Appendix |A] and in the caption. The decrease in pq is mostly due to the increase 
of t°(x) = T°(x)/T°(0) (dash-dotted curve in Fig. [SJ. On the other hand pu has to increase with x in order to 
achieve the crossover from T c enhancement (T c (0)/T°(0)=2.5) at x=0 to T c depression (T c (l)/T°(l)=0.34) at 
x=l. This observed crossover from pair forming to pair breaking by CEF excitations is only possible if the dipolar 
scattering strength increases with x. This is in part due to the increase of the tetrahedral CEF parameter ycF 
(Appendix lAl since for small ycF the dipolar matrix element d ~ y^ F shows a strong increase with ycF- The 
latter is compensated by the t°(x) increase leading to an almost flat pd(x) for large x. The associated crossover 



10 



from T c (a;)/T° enhancement to reduction obtained from Eq. (l36l) is shown in Fig. [5] with an almost fiat reduction 
factor at large x. Altogether, because the background T®(x) strongly increases with x, the renormalised T c (x) 
exhibits a minimum for intermediate concentrations (bottom Fig. [5j. The experimental values from Ref. are 
shown for comparison (squares). The T c (x) curve was calculated for 7 = and moderate 7 = 0.15 with only small 
difference, especially for larger x (around the minimum region). The precise form of Pq{x) and pd{%) cannot be 
determined presently because no reliable information on CEF parameters and matrix elements for intermediate 
x is available. However the crossover from mainly pair forming at small x to pair breaking behaviour at large x 
(FigU) is robust. We conclude that the vibronic splitting does not play an essential role in the T c (x) minimum 
formation. This is due to the fact that for x ~ x c close to the crossing region already Sl s (x c ) /2T C 3> 1 where the 
pair breaking functions Sx 2(2^) (Appendix IB"]) vary slowly with x s . Therefore the coupling 7 > hardly affects T c 
for larger x. Its effect would be much bigger if the mode crossing would appear at energies comparable to T c , i.e. 
il s (x)/2T c ~ 1. In the case of Pr(Osi_ x Ru x )4Sbi2 this is not possible because already for x=0 we have Ao/2T° = 5.4. 

One may conclude that the T c (x) minimum is not directly linked to the crossing of rattling phonon mode and 
CEF excitation found in Ref. [2(J It is rather a combined effect involving the crossover from T c /T° enhancement 
to reduction and and the increase in the background T^(x). In this scenario the observed T c (x) minimum also does 
not imply or suggest a symmetry change of the order parameter from below to above the Ru concentration at the 
minimum. 



VI. SUMMARY AND CONCLUSION 

In this work the possible nonadiabatic effects of CEF singlet-triplet excitations and rattling phonons of rare earth 
hosts in the ca ges of Pr(Osi_ x Ru x )4Sbi2 have been investigated. This is suggested by specific heat experiments of 
Miyazaki et al [201 ] which show a crossing of triplet excitation and rattling phonon energies at an intermediate Ru 
content. 

It has been proposed that a magnetoelastic coupling between the singlet-triplet excitations and the local rattling 
modes should lead to vibronic splitting and mixed-mode formation around the crossing point. These features can be 
detected in the spectral functions measured by IKS experiments. It should also be observable in the low temperature 
depression of the symmetry elastic constant as function of Ru concentration which measures directly the quadrupolar 
susceptibility of the vibronic excitations. On the other hand the magnetic susceptibility is not affected by the mode 
splitting. Its comparison with experiment indicates an increase of dipolar matrix elements for increasing x and 
possibly the influence of the higher lying triplet r^. 

The superconducting T c (x) behaviour of Pr(Osi_ a; Ru K )4Sbi2 shows an enhancement at small and reduction at 
larger Ru concentration compared to the background T°(x) of La(Osi_xRu x )4Sbi2. It also exhibits a minimum 
at intermediate concentration. The analysis presented here suggest that this behaviour is the result of a crossover 
between primarily quadrupolar pair formation and mainly dipolar pair breaking mechanism originating in the singlet- 
triplet excitations. The vibronic coupling has little influence on the existence of the minimum (Fig. [5]) because of 
the large singlet-triplet splitting for x ~ 0.65 in comparison to T c . The minimum is rather a combined effect of a 
decreasing enhancement factor and the increasing background T°(x). An essential feature of the model is a growing 
dipolar (magnetic) character of the triplet for larger Ru content which may be due to an increase of the tetrahedral 
crystal field part. For a better determination of the model parameters it would therefore be important to determine 
the CEF parameters for intermediate Ru content by INS experiments. The experimental knowledge of T®(x) of 
La(Osi_ a; Ru x )4Sbi2 in the whole Ru concentration range would also allow further improvement of the model. 

Appendix A 

In this appendix we summarize basic properties of the tetrahedral CEF model and interaction model needed in 
the analysis. The CEF potential in tetrahedral symmetry is given in Ref. Hi]. Aside from an overall scale W it 
is determined by two parameters xcf(%), Ycf(x) which characterize fourfold cubic and tetrahedral contributions 
respectively. These parameters will depend on the Ru content x. The same is then true for dipolar and quadrupolar 
matrix elements b 2 D (x) = (20/3)d 2 (x) and b 2 Q (x) = (140/3)(1 - d 2 (x)) which depend on a mixing parameter d given 
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by [22 



{X> 2 I (3 + 2xcF(x))2 + 1008ycF(xr) [ ' 

.(2) 



It is a measure of the tetrahedral CEF since d — for ycF = 0. In this model the singlet-triplet ri — T 4 splitting is 
given by 22] 

A(x) = 2VF(36 - 58x C f(x)) - 4VF[(3 + 2x C f{x)) 2 + 1008y£ F (x)] * (A2) 

For PrOs 4 Sbi2 (x=0) we use x C f(0)=0.45 and j/cf(Q) = 0.1 [22]. In PrRu 4 Sbi 2 (x=l) the CEF splitting is much 
larger. This suggests that xcf(1) is close to zero according to the LLW tables [3a]. Furthermore for x = 1 the dipolar 
matrix element b 2 D and hence ycF has to increase. Therefore we use the set xcf(1) = 0.0 and J/cf(1) = 0.2. It leads 

(21 (2) 

to comparable intensities for the ri — > T\' (84 K) and ri — > T\' (145 K) transitions in qualitative agreement with 
INS results in Ref. [37|. Furthermore the CEF scale factor in Eq.([A"2"]) is given by W ~ 1.9 K leading to A(0) = 8K 
and A(l) = 84 K. For general x an interpolation between values at x = 0, 1 is employed according to 

Xcf{x) = XQpil — x) + Xq F x + Xq F x(1 — x) (A3) 

and similar for ycF(x). The interpolation parameter sets (xq F , Xq F , x 2 cf ) = (0.45,0,0) and {yc F , J/cf> Vcf) = 
(0.1,0.2,0.1) have been used for calculating d 2 (x) from Eq. (|A1[) and then b 2 (x)) (i — Q,D) needed in Eq. (|37l) . 
These CEF parameters also reproduce the experimental A(x) behaviour. 

For the calculation of pair forming and breaking functions in Eq. (|37[) one needs the dimensionless interaction 
constants pi(x) (i = Q, D) in additon to the CEF matrix elements. The former are obtained by a similar interpolation 

as in Eq. (TM| with parameter sets (p°q,Pq,Pq) = (22.34,22.34,0) and (p° Q , p x Q , p 2 Q ) = (3.2,10.3,6.0). The p°' 1} are 
chosen such that for x = 0, 1 the absolute value of T c (x) corresponds to the experimental one in Fig. [SJ Finally the 
p 2 are adjusted to obtain the minimum position and value of T c (x). 

Appendix B 

The functions Si^{x) in Eq. (13"51) were derived in Ref. [33| and are given here for completeness: 

where ip{ z ) = d\nT(z)/dz is the digamma function [36| 
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